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splitting of charge carriers in solids. While many of our results apply more generally, we focus 
especially on the Rashba model describing electrons in the conduction band of asymmetric 
semiconductor heterostructures. Appropriate semiclassical limits are distinguished that 
describe various situations of experimental interest. Our analytical fomulae, which take full 
account of Zeeman splitting, are used to analyse recent magnetic-focusing data. Surprisingly, 
it turns out that the Rashba effect can dominate the splitting of cyclotron orbits even when the 
Rashba and Zeeman spin-splitting energies are of the same order. We also find that the origin 
of spin-dependent cyclotron motion can be traced back to Zitterbewegung-like oscillatory 
dynamics of charge carriers from spin-split bands. The relation between the two phenomena 
is discussed, and we estimate the effect of Zitterbewegung-related corrections to the charge 
carriers' canonical position. 
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1. Introduction 

Magnetic focusing of ballistic charge carriers in metals [1,2] and semiconductors [3-5] has 
been successfully used to elucidate fundamental materials properties such as the shape of the 
Fermi surface [6, 7], Andreev reflection in superconductor/normal-metal hybrid structures [8, 
9], surface crystallography [10], phase coherent transport [3], and emergent quasiparticles 
in the fractional-quantum-Hall regime [11]. The fundamental setup of a magnetic-focusing 
experiment is quite simple; see Figure 1 . It requires sufficiently ballistic transport between 
two fixed, co-linear (injector and collector) contacts. A large number of injected-particle 
trajectories converge at the collector every time the contact separation L equals an integer 
multiple of the cyclotron diameter 2rc. Peaks occurring in the measured collector voltage 
at concomitant magnetic fields are the experimental signature for magnetic focusing. 
Recently, magnetic-focusing trajectories of electrons in a two-dimensional semiconductor 
heterostructure have been imaged directly using scanning-probe microscopy [12]. 

Several experiments [13-17] have investigated the possibility to spatially separate charge 
carriers belonging to spin-split bands using the magnetic-focusing technique, which may have 
ramifications for the emerging field of spin electronics [18, 19]. The desire to use magnetic- 
field-independent, spin-orbit-induced spin splitting for this purpose [14-17] is fuelling 
renewed interest [20-24] in the theoretical study of spin-orbit effects in the semiclassical 
regime [25-37]. We present a critical analysis of the interplay between spin-orbit coupling 
and cyclotron motion in both quantum and semiclassical regimes, with particular emphasis 
on spin-dependent magnetic focusing. We also discuss this effect in the context of another 
topic of great current interest, namely Zitterbewegung in solid-state systems [38-44], and 
comment on the (un-) suitability of proposed intuitive interpretations in terms of a spin- 
dependent focusing field [14] or a spin-dependent Lorentz force [45-47]. 

This article is organised as follows. We start by reviewing semiclassical theories that 



Figure 1. Schematic setup for a magnetic -focusing experiment. Current is passed through 
the injector contact, and the voltage in the collector contact is monitored as a function of a 
magnetic field applied in the z direction. The latter forces charge carriers to move on cyclotron 
orbits in the xy plane. Particles injected on trajectories starting out sufficiently close to parallel 
to the y direction will be focused into the collector contact when the contact separation equals 
an integer multiple of their cyclotron-orbit diameter. As a result, peaks are observed in the 
collector voltage at the corresponding magnetic-field values. 
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have been developed for systems with finite spin spUtting. These theoretical approaches 
are then applied to describe cyclotron motion of charge carriers in two-dimensional (2D) 
heterostructures subject to both Zeeman and Rashba [48, 49] spin splittings. Our results are 
used to analyse recent experimental data obtained from p-type GaAs [14] and n-type InSb [16, 
17] samples. Subsequently, we discuss the intricate connection between spin-dependent 
magnetic focusing and Zitterbewegung of charge carriers from spin-split bands [39-41,43]. 
Conclusions are given in the final Section. 

2. Basic aspects of semiclassics in presence of spin splitting 

Non-relativistic single-particle Hamiltonians with spin splitting can generally be written in 
the form [32-36] 



Here we denote the particle's position, momentum, and spin angular-momentum operators 
by r, p, and S = /Is, respectively, where s is the vector of generators for SU(2) rotations 
in the spin-s representation, and h is the Planck constant. The vector operator B represents 
an effective magnetic field that contains, in general, contributions due to the Zeeman effect 
and spin-orbit coupling. Note that B has the dimensionality 1/time. In the following, we 
distinguish three possible semiclassical limits. To keep our notation uncluttered, we do 
not explicitly distinguish quantum-mechanical operators from their associated semiclassical 
phase-space symbols. We will clearly separate the discussion of purely quantum and 
semiclassical properties to avoid any possible confusion arising from this simplification. 

2.1. Precessing-spin semiclassics 

Representing a truly quantum-mechanical correction, the second (spin-splitting) term in (1) 
vanishes in the usual semiclassical limit that is defined as — > such that | s | = \S\/h remains 
constant. As a result, the semiclassical orbital dynamics is unaffected by the spin degree of 
freedom [31-33]. The classical trajectory {r{t), p(i)}, as determined by Hq, prescribes the 
dynamics of the classical spin s via a precession-type equation of motion 



i.e., we can interpret \B\ as the precession frequency. This type of semiclassics has also been 
called a weak-coupling limit [32-36]. Physically, it corresponds to the regime of a perfectly 
classical, spin-independent orbital motion with an associated trajectory-dependent precession 
of a classical spin [31-34]. 

In many experimentally relevant situations the appropriate semiclassical description 
will be of the precessional type discussed here. This has, e.g., been shown for anomalous 
magneto-oscillations [37]. However, there are experimentally accessible regimes where 
the spin dynamics actually influences the orbital motion. Spin-dependent cyclotron motion 
represents a pertinent example [14, 16, 17]. The semiclassical description of such situations 
is desirable, motivating the consideration of alternative schemes for performing the classical 



7i = 7io(r,p) + B(r,p)-S 



(1) 



s = B(r(t),p(i)) xs 



(2) 
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limit. We proceed to discuss two of these, both of which correspond to a strong-coupUng-type 
semiclassics. 

2.2. Spin-orbit-intertwined semiclassics 

One possibility to keep the spin-splitting term in the Hamiltonian (1) finite in the semiclassical 
limit — i> is to simultaneously require |S| to remain constant, which implies |s| — > oo. The 
time evolution of the spin state will then affect orbital dynamics and vice versa. The resulting 
set of semiclassical equations of motion is given by [33] 



Hence, in this case, the spin and orbital dynamics are mutually affecting each other. These 
equations of motion have previously occured in [35, 36], but without noticing that they only 
provide the leading semiclassical dynamics when |s| ^ cxo is considered in addition to h ^ 0. 
Moreover, in [35, 36] the notion of an extended phase space was introduced, intending to 
stress the role of spin as an independent classical dynamical variable. Kinematically, the same 
phase space arises in the precessing-spin semiclassics; however, in that context the spin-orbit 
dynamics are not of a Hamiltonian form because the spin dynamics is driven by the orbital 
motion, without any feedback of the spin dynamics on the orbital motion. 

2.3. Adiabatic-spin semiclassics 

An alternative, in some sense very-strong-coupling [25-27,32,34,36] semiclassics is 
obtained by letting /i — > while keeping the spin projection frozen to a quantised value 
Sz = with respect to the instantaneous direction of 3. Then the associated orbital 
dynamics is governed by the Hamiltonian 



For each possible value S^, ranging over S, S — h, . . . , —S, a generally different classical 
trajectory is obtained. In contrast to the two cases discussed above, here the spin degree of 
freedom itself has no dynamics; it just introduces a Berry-phase-like contribution to the orbital 
motion [50]. To avoid the problem of mode-conversion, this approach is restricted to the case 
|B(r, p) I > for every point along a trajectory [25, 26]. 

3. Application to the Landau-Rashba model 

To be specific, we consider the Landau-Rashba model [48,49] that describes spin s = 1/2 
conduction-band electrons in an asymmetric 2D heterostructure subject to a perpendicular 
magnetic field B = Vr x A = z. It is of the form given in (1), with 



P = - VrU , 

S = B(r,p) X S 



(3a) 
(3b) 
(3c) 



ns^^no{r,p) + S,\B{r,p) 



(4) 




(5a) 



(5b) 
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Here m and g are the effective mass and Lande factor of the 2D electrons, c^co = eB/niQ 
with mo the electron mass in vacuum, and a characterises the strength of the Rashba spin 
splitting [48,49]. The 2D hetero structure growth direction is taken as the Cartesian z axis; 
with z being the associated unit vector. 



3.1. Quantum solution: Jaynes-Cummings model 

A complete quantum solution of the Landau-Rashba model is available [51], as it is 
equivalent [36] to the exactly soluble Jaynes-Cummings (JC) model [52] in the rotating- 
wave approximation. (The JC model describes coupling of a harmonic oscillator, here 
associated with the Landau levels, to a two-level system, here represented by the electronic 
spin degree of freedom.) Using techniques developed in theoretical quantum optics [53], we 
recently obtained [43] the exact Heisenberg time evolution of spin and position operators. 
The perpendicular-to-the-plane spin component can be separated into two parts, Sz{t) — 
Sz + Sz{t), with time-independent and oscillating parts given by 

' ~4BJC-S~2l2 my2BJ°-S ' ^ ' 

Sz{t) = {Sz - S,) exp {-2it • S/h) . (6b) 

Here = B(r,p) — {mo / m)uJco z is an effective magnetic-field operator that governs 
spin precession in the Landau-Rashba model. Similarly, the 2D position operator can be 
decomposed into a constant part, which corresponds to the guiding-centre position of a 
cyclotron orbit, and a time-dependent oscillatory part. For the sake of brevity, we will 
use a compact complex notation [54] for 2D position r = (x, y) and kinetic-momentum 
TT = p + eA = (ttx-, Tiy): R = x — iy,Ii = 'Kx — iny, and extend this notation also to the 
in-plane spin components: S± = (Sx ± iSy) /2. In the Heisenberg picture, the time evolution 
of the complex 2D position is then given by R{t) = R + R{t), with 

in 

R , (7a) 



JC 



R{t) — exp 



-it u, 



, , iU smicust) in B^^ ■ S 
cosfa;^^) \ : — 

m 

Here ujs = \J (B''^ ■S/Kf' + hjjJcma^/2, and ujc = eB/m is the cyclotron frequency of 
electrons in the semiconductor material. 

In principle, the expressions given in (6a)-(lb) allow for the calculation of time- 
dependent spin and position expectation values for any initial state. In practice, such a 
calculation may turn out to be rather cumbersome and difficult to interpret. Hence, in the 
following, we will discuss magnetic focusing in the context of semiclassical approaches 
applied to the Landau-Rashba model. Approximate semiclassical approaches are often 
practical for understanding certain physical phenomena and also provide a rather general 
framework to treat quantum systems of interest. However, the above exact results provide 
a useful benchmark for their reliability. 
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3.2. Precessing-spin semiclassics of the Landau-Rashba model 

In this case, the orbital dynamics is entirely governed by the Landau model. Using the 
compact complex notation introduced above, we have 

Rit) = R+ 6-^"^"=* , (8a) 
miUc 

n(i) =ne-^'"'=* . (Sb) 

The spin dynamics is determined by the precession equation (2), where B is given by (5b). 
For the z and in-plane components and s± — (s^ ± isy) /2, they read explicitly 

s^^ --U{t)sz{t)-i^u,os4t) , (9a) 

=a[s4t)U*{t) + s+{t)U{t)] . (9b) 

[The equation for follows from complex conjugation of (9a).] From (9a) and (9b), it 
follows that 

Sz = —ujp Sz , (10a) 
with the spin-precession-related frequency scale 



u;p = ^/(a7r)^ + a;2(^l-^j . (lOb) 

Note that TT = |7r| = |n| = y/2mE. Integration yields 

Sz = C+ e^'^P* + C_ e-^'^p* , (11a) 

s,{t) = Sz-C+ — + C_ : , (lib) 

lUJp tOJp 

with arbitrary constants C±. Inserting this result into (9a) enables one to find S-(t). We omit 
this step here. 

It is illustrating to note that hujp/2 emerges as the eigenvalue of the JC Hamiltonian 
B''*" ■ S in the limit where the kinetic-momentum operators tt^ and Tiy are treated as c- 
numbers. This corresponds to the early semiclassical treatments of the JC model [55]. Thus 
(1 IZ?) with (lOb) reflects the exact quantum solution for spin precession in the Landau-Rashba 
model (6b) taken in the appropriate weak-coupling limit. Our results obtained here generalise 
those presented in Sec. 5.1.1 of Ref. [36] where the weak-coupling limit of the Landau- 
Rashba model was previously discussed. The complete disappearance of h from the dynamics 
described in this Section is an expected feature of the precessing-spin semiclassical limit. 

3.3. Adiabatic-spin semiclassics of the Landau-Rashba model 

In the limit where the electron spin is assumed to be either aligned or anti-aligned with the 

local field B(r, p) [given by (5b)], the Landau-Rashba model specialises to a pair of terms of 
the form (4) that govern the dynamics of electrons with spin projection ±1/2. As the Rashba 
and Zeeman contributions to B(r, p) are perpendicular to each other, we have 



B(r,p)| = J(|a;eo) + aH^ . (12) 
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Thus TT = |7r| = \/2mE is again a constant of the motion. For a fixed value of conserved 
energy E, it assumes two different values 7r<j for particles distinguished by spin projection 
cr/2, where a — ±1. The values of tt^ can be found from 



The equations of motion resulting from Hamiltonians (4) for the Landau-Rashba case 
can be written as 

r = — , (14a) 
-n- = - er X B , {Ub) 

which describe cyclotron motion with a spin-dependent effective mass 

rucr . mha/2 

— = 1 - (7 , , (14c) 

^ ^J2mE + {mhajiy + {gu^^l{2a)f 

and thus spin-dependent frequency ujccr — eB/m^. The cyclotron radius Tc = T^/{muJc) is 
different for the two spin species because of their different values of tt = tTo^ for fixed energy 
E. A straightforward calculation yields 



^ccr 

mujr 



2mE + 2 ' - a haJ2mE + (m^) ' + f ^) '.(15) 

iLJc \ \2/ y \2/V2q;/ 



A few comments about our results are in order. Firstly, the fact that the above expressions 
for spin-dependent cyclotron frequency and radius depend on the parameter ha/ 2 (which has 
the dimension of velocity) is a direct consequence of the way the adiabatic-spin semiclassical 
limit is performed. Secondly, for = 0, our expression for tUca is exactly the same as 
that found in Ref. [21] where, ostensibly, the spin-orbit-intertwined semiclassical limit was 
discussed. Also, our result for agrees with the corresponding expression from Ref. [21] 
to leading order in the large-£' limit. Apparently, the approximate scheme employed by the 
authors of Ref. [21] is essentially equivalent to the adiabatic-spin semiclassics, even though 
they recover a finite z component of spin as given in Equation (6a) above. Thus a consistent 
treatment of the Landau-Rashba model using spin-orbit-intertwined semiclassics seems to be 
still lacking. Thirdly, the canonical equations of motion (I4a) and (\4b) indicate that the effect 
of adiabatic-spin semiclassics on the orbital dynamics is better described as a renormalisation 
of the effective mass [22,48] than a renormalisation of the focusing field [14]. Lastly, (14^?) 
represents the familiar expression of the Lorentz force in terms of a particle's velocity without 
any trace of the previously claimed [46] spin-dependent contribution. Such a result is expected 
from proper quantum-mechanical derivations [56] of the Lorentz-force operator. 



3.4. Analysis of magnetic-focusing experiments 



To leading order, no spin splitting of magnetic-focusing peaks occurs in the precessing-spin 
semiclassical limit. For this case, back action of spin dynamics on orbital motion may appear 
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Table 1. Parameters associated with and extracted from recent spin-dependent magnetic- 
focusing experiments. Besides contact separation L and focusing fields B±, we also provide 
the value fcpo = \/27rno of Fermi wave vector as derived from the 2D sheet density no, which 
can be compared to kp extracted from the focusing data. 



material 


gm 

2mo 


L [nm] 


fcpo [nm 1] 


B+ [T] 


B- [T] 


fcp [nm 


fcso [nm 1] 


p-GaAs^ 


1.4 


800 


0.093 


0.17 


0.20 


0.11 


0.009 


n-InSb* 


0.36 


600 


0.14 


0.30 


0.36 


0.15 


0.014 


n-InSb* 


0.36 


600 


0.19 


0.37 


0.50 


0.20 


0.030 



^ From Ref. [14]. As in this work, we apply the fc-Unear Rashba model for conduction- 
band electrons to interpret the data. A more detailed study would have to take into account 
fundamental differences between Rashba spin sphtting in 2D electron and hole systems [58]. 
* From Ref. [16]. 



* From Ref. [17]. The given value of A;po is derived from the focusing field (0.42 T) expected 
in the absence of spin splitting (as stated by the authors). 



only in corrections of order h that can be included in principle [57]. Instead of considering 
this possibility, we focus here on the adiabatic-spin semiclassics which provides a proper 
description of magnetic focusing already in leading order. 

It is useful to define effective wave- vector scales k-p = ^/2mE + {mha/2)'^/h and 
^so — iTia/2, which are associated with the 2D carrier sheet density and Rashba spin splitting, 
respectively. After equating 2rca from (15) with the contact separation L and performing 
some straightforward algebra, we find a relation that has to be satisfied by each of the two 
experimentally observed focusing fields B±: 

+ j = [^^J .(16) 

Thus measurement of the spin-split focusing peaks allows to determine both kp and kso 
directly, with the effective electron mass entering only the parameter gm/ (2too). For g — 0, 
the relation (16) specialises to eLB^j = 2h{kF — akso), which was used in Ref. [14] to extract 
kso in a GaAs 2D hole system. Table 1 summarises values obtained from existing magnetic- 
focusing data, taking into account the finite Zeeman splitting. It turns out that, for the sample 
parameters realised in these experiments, the dependence of extracted k^o on gm/{2mo) is 
rather weak over an extended range before Zeeman splitting becomes suddenly dominant. 
This surprising feature, which is illustrated in Figure 2, explains why it was possible to extract 
a reasonable value for kso in Ref. [14] even though, in that experiment, Rashba and Zeeman 
spin splittings were of comparable magnitude for states at the Fermi energy. Note, however, 
that the range of parameter gm/ (2mo) over which Rashba splitting can be reliably extracted 
will be reduced in samples with smaller contact separation and concomitantly higher focusing 
fields. 

The values of kgo given in Table 1 are on the order of 10% of the effective Fermi wave 
number kp. Thus the applicability of the adiabatic-spin semiclassics for typical experimental 
situations could be questioned. A detailed discussion of this point would benefit from a fuller 
understanding of spin-orbit-intertwined semiclassics in the Landau-Rashba model, which is 
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gm 
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Figure 2. Dependence of extracted kso on the assumed value for the reduced g-factor 
gm/{2mo). The calculation of the solid red (dashed blue) curve used focusing data from 
Ref. [16] (Ref. [14]). Apparently, a rather weak variation of extracted fcgo with gm/{2mo) 
persists to quite large values of the latter which, in experiment, correspond to comparable 
magnitudes of Rashba and Zeeman spin splittings. 



currently lacking. Incidentally, results from a numerical simulation [20] performed using 
experimentally realistic parameters provide strong support for the assumption of adiabatic- 
spin dynamics. 

Our analysis of experimental data focused exclusively on the first magnetic-focusing 
peak. This allowed us to neglect scattering at the lithographic barrier between injector and 
collector contacts, which is relevant for higher-order focusing peaks. Within a ballistic 
semiclassical approach such as ours, spin flips occuring during collision with the sample edge 
can be taken into account phenomenologically [59] by including the possibility for particles 
to continue on either one of the spin-split cyclotron orbits after each reflection. This model 
predicts that the second peak will be unsplit (split into three parts) in the absence (presence) 
of boundary spin-flip scattering [17]. 

4. Relation to Zitterbewegung 

Zitterbewegung (ZB) was originally introduced by Schrodinger as the technical term for an 
oscillatory orbital motion performed by free relativistic electrons whose dynamics is governed 
by the Dirac equation. See Refs. [60-63] for modern descriptions of the effect. ZB has never 
been directly observed, partly because the associated period and amplitude (< 10^^^ s and 
< 0.004 A, respectively, for electrons in vacuum) are out of reach for any current experimental 
equipment. For ultra-relativistic particles, the ZB amplitude is of the order of the de Broglie 
wave length [43], limiting the suitability of scattering experiments to detect the effect. 

Analogs of ZB in a nonrelativistic solid-state context have recently attracted great 
interest [38-44]. In particular, charge carriers from spin-split bands are expected to perform 
an oscillatory motion that is entirely analogous to ZB [39-41,43]. The experimentally 
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observed [64, 65] zero-field spin precession of electrons and holes turns out to be closely 
related to ZB [43], but no direct ramification of ZB in coordinate space has been measured. 
Theoretical studies suggest [41,43,66] that ZB results in a spatial separation of carriers 
with opposite spin that is of the order of the de Broglie wave length. Here we show that 
spin-dependent magnetic focusing is also closely related with ZB. For the sake of notational 
simplicity, we neglect Zeeman splitting from now on and consider only the limit of sufficiently 
small magnetic fields where Landau-quantisation effects are negligible. 



4.1. Cyclotron motion of charge carriers performing Zitterbewegung 

It is a consequence of ZB that the time-dependent position [velocity, spin] operators r (t) [v(t) , 
S(i)] can be written as the sum of an average (smoothened over ZB) part f(t) [v(t), S{t)] and 
an oscillatory part r{t) [v{t), S{t)]. Universal expressions, in terms of suitably defined ZB 
frequency and amplitude operators a)(p) and F, have been obtained for all these operators for 
a range of (multi-band) models [43]. For example, the average part of the velocity is given by 

.^f^-F (17) 
dp 

for any two-band Hamiltonian H, including the Rashba model. We continue by discussing 
this special case only. A finite acceleration of the (ostensibly free!) particles performing ZB 
in zero magnetic field is found 

VzBW=i^(p)Fe-^'^(P)* = a25,(i) pxz , (18) 

where the r.h.s expression is the specialisation to the Rashba-model case, and we used the 
appropriate expression for F that, within our current notation, reads 



d{B-S) fhaV p 



dp V 2 y B • S ■ ^^^^ 
As the r.h.s of (18) indicates, this acceleration is intimately related to spin precession in the 
Rashba model. In particular, the fact that S*^ = in zero magnetic field is directly associated 
with the vanishing time average of the ZB-related acceleration (18) for this case. 

A finite perpendicular magnetic field forces charged particles on cyclotron orbits and, 
therefore, leads to additional time dependences of their dynamical variables. To elucidate 
the interplay between ZB and classical cyclotron motion more clearly than it emerges, eg, 
from the exact solution of the JC model given in Sec. 3.1, we concentrate on the low-field 
limit where the time scales associated with these two effects are well-separated. This regime 
allows one to consider quantities that are averaged over the ZB time scale but are still time- 
dependent because of the cyclotron motion. 

As indicated by (6a), Sz becomes finite in a perpendicular magnetic field. In that 
situation, the average of the acceleration (18) over the ZB time scale becomes finite and 
yields, in the low-field regime and with Zeeman splitting neglected, the Lorentz-force-like 
term 

VzB(i) = -- U{t) - ^) X B . (20) 
m \ m 
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To obtain (20), we used {6a) above, as well as the relation (17) specialised to the Landau- 
Rashba model in the low-field limit (where the non-commutativity of t^^ and tTj, can be 
neglected), 

- + ^U . (21) 

The total ZB-averaged acceleration experienced by a charged particle subject to a magnetic 
field is given, in the low-field limit, by the sum of the ordinary Lorentz-force contribution and 
the finite ZB-averaged acceleration (20). It turns out to have the form 

if{t) = -—Y{t)xB , (22) 
m 

ie, is determined by the ZB-averaged velocity v{t) derived from (21). Applying the adiabatic- 
spin semiclassical limit, (21) specialises to (14fl), and (22) is equivalent to (14-b). Hence 
the spin-dependent cyclotron mass m^^ emerges because of the ZB contribution (20) to the 
total acceleration (22). Similar to the spatial separation of spin-polarised partial waves for an 
unpolarised electron beam injected into a wave guide [66], spin-dependent cyclotron motion 
is thus a direct consequence of ZB -related dynamics arising in the presence of spin splitting. 

As noted already a long time ago [61], ZB -related terms can contribute to expectation 
values of observables such as (r^). Similarly, the evaluation of (r^) for eigenstates of the 
Landau-Rashba model at fixed energy E yields a hint of the existence of two different 
cyclotron orbits for electrons from spin-split bands [21]. 



4.2. Corrections due to anomalous position operator 

The effective Rashba Hamiltonian describing electrons in the conduction band can be thought 
of as arising from a canonical transformation (Lowdin partitioning [58]) that is similar to the 
Foldy-Wouthuysen transformation [67] needed to arrive at the proper non-relativistic limit of 
Dirac-electron theory. That same transformation will change the form of the physical position 
operator associated with a point-like particle, which will then be different from the canonical 
position operator acting in the reduced Hilbert space of the conduction band. For Dirac 
electrons, the two position operators have been proposed to be associated with a particle's 
centre of charge and centre of mass, respectively [61]. See Ref. [68] for a closely related 
discussion refering to solid state systems. Transport measurements such as the ones employed 
by the magnetic-focusing technique can be expected to be sensitive to the physical (charge) 
position rather than the canonical (mass) one. 

The formal relation between physical position Tp and the canonical one r is 

rp = r + — TTxS , (23) 

where A is the (effective) Compton wave length. Typical A-values for carriers in generic 
semiconductors can be found in Refs. [38,43]. In a finite magnetic field, the shift between 
physical and canonical position gives rise to a correction to spin-split cyclotron radii that is of 
the order of {eBlA^/h and, therefore, usually quite small. A possible exception could be InSb 
where A 4 nm. 
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5. Conclusions 

We have studied theoretically the cyclotron motion of charge carriers from spin-split 
bands, highlighting exact quantum and semiclassical results for electrons in asymmetric 2D 
semiconductor heterostructures. A spin-dependent splitting of cyclotron orbits is found in 
the adiabatic-spin semiclassical limit. Relevant parameters of real samples used in recent 
magnetic-focusing experiments were extracted, taking full account of Zeeman splitting. Our 
analytical formulae should also be useful for analysis and design of future spin-dependent 
focusing measurements. We furthermore elucidated the intricate relationship between spin- 
split cyclotron orbits and Zitterbewegung of charge carriers in systems with strong spin-orbit 
coupling. 

Future studies will be aimed at a systematic investigation of similarities and differences 
exhibited in the cyclotron motion of particles from generic two-band models such as those 
describing relativistic Dirac electrons [69] or holes in a typical semiconductor's valence 
band [70]. 
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